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Background 

Cardiac trabeculae are complex and difficult to quantify 
but measuring their biological signal in many cardiac 
diseases (left ventricular [LV] non-compaction, hyper- 
trophic cardiomyopathy [HCM]) has value. We have 
recently shown how trabeculae differ between ethnicities 
in health and between HCM mutation carriers without 
hypertrophy and healthy volunteers. Here we present a 
novel tool for measurement for LV trabecular complex- 
ity by cardiovascular magnetic resonance (CMR) using a 
fractal analysis. Theoretical fundamentals are discussed. 
They include the measurement procedure, mathematical 
modeling and performance evaluation through validation 
by digital phantoms and reproducibility analysis when 
applied to real-world cardiac data. 

Methods 

The measurement system consists of an in-house devel- 
oped software written in MALTAB R2012b that para- 
metrizes 2-dimensional LV cine end-diastolic short-axis 
frames in DICOM image format. Images are scaled, 
masked and segmented before undergoing fractal analy- 
sis. Endocardial contours serve as the basis for estima- 
tion of the fractal dimension (FD) -a unitless index of 
endocardial complexity, calculated using a box-counting 
algorithm which cycles through a series of grids of 
boxes overlying the region of interest. Validation was by 
digital phantoms consisting of planar targets of known 
FD (Figure 1) generated in MATLAB. The procedure 
was further evaluated in terms of intra- and inter-obser- 
ver variability of calculated maximal apical FD in 45 
randomly-selected scans from within a research CMR 



database (all 1.5-T, Siemens and General Electric) and 
in terms of inter-study reproducibility (20 repeat scans, 
same 1.5T magnet; mixed healthy/diseased cohort). 
Impact of varying slice thickness (43 cine slices repeated 
at thicknesses: 8, 7, 6 mm) and magnetic field strength 
(134 cine slices repeated at 1.5 and 3T; mixed healthy/ 
diseased cohort) on raw slice-by-slice FD was also 
evaluated. 

Results 

Phantom validation of the box-counting algorithm pro- 
duced measured FD which closely approximated 
expected values as proof of reliability. Intra- and inter- 
observer variability for maximal apical FD was high 
(ICCs: 0.94, 95%CI 0.90-0.97 and 0.92, 95%CI 0.86-0.95 
respectively; Figure 2A, B) as was inter-study reproduci- 
bility (ICC 0.96, 95%CI 0.83-0.98; Figure 2C). ANOVA 
with repeated measures determined that mean FD did 
not differ significantly between the 3 slice thicknesses 
(F(2, 84) = 2.259, p = 0.1). Agreement between raw FD 
values was good irrespective of magnetic field strength 
(ICC 0.92, 95%CI 0.90 0.94; Figure 2D). 

Conclusions 

This is the first comprehensive implementation of a 
fractal procedure applied to the clinical CMR domain. 
Fractal analysis provides a way of collapsing the mor- 
phometric complexity of LV trabeculae into a single 
numerical value which can be compared between groups 
of subjects. Here we show how fractal analysis is a 
robust, reliable and reproducible method for capturing 
the richness of cardiac trabeculae imaged by CMR. 
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Cardiac Contour with Meshgrid Overlay in MATLAB 



O OOOO:0 

o oooo o 
o 0000:0 
o 0000:0 
o 0000:0 
o 0000 0 
o 0000:0 
O O00O: 0 

" — ^00:0 




0000; 
0000: 
0000 
0000 



, 0000:0 
'o 0000: 0 
o 0000:0 
o 0000:0 
o 0000:0 
o 0000:0 



OOOO:' 
oooo:' 
oooo:' 
0000:' 

OOOO' 

0000:' 

OOOO'' 
OOOO:' 

0000:' 
OOOO': 
OOOO :' 
OOOO ' 
OOOO': 
OOOO : 

OOOOj 

0000:' 
OOOO :> 
0000:' 

OOOO'l' 

0000:' 

OOOO'-' 
OOOO ' 
OOOO :' 
OOOO : 

OOOO' 

0000:' 
0000:' 
0000:' 
0000:' 
OOOO :■ 



Figure 1 Using the Pascal analogy triangle, we constructed mathematical fractal phantoms: Sierpinski gaskets of px heights:! 6-32-64- 
128-256 and known FD. 15, 32 and 64 px gaskets are reproduced here(left). Predicted/measured FD demonstrate high concordance. To the 
right is a cardiac contour after extraction from DICOM data In MATLAB, undergoing fractal analysis. Each slice is scaled using a magnification 
algorithm with bicubic interpolation{Stanescu E et al.), masked and then segmented using a region-based level set method(Li C. et al). The 
contour, composed of N pairs of positive long X/Y coordinates, is overlaid by a series of meshgrids, replicated over 4 random start points. 
Starting box-size = 45% that of the minimum circumscribed box with parallel constraint; terminal box-size = 4 px to maintain stable regression. 
Box size decreases by one px per iteration. At every level, scale and box count are recorded and formulated into a natural logarithmic plot. The 
slope determines the FD. Px = pixel. 
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Figure 2 Scatter plots showing difference in FD measurements against average of the repeated values for: A) intra-observer 
variability; B) inter-observer variability; C) inter-study reproducibility; D) inter-field reproducibility. There is no systematic bias Middle 
ine = mean difference, two extreme lines = limits of agreement (+1.96 and -1.96 standard deviation). 
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